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' We prove that static, spherically symmetric, asymptotically flat, 

regular solutions of the Einstein- Yang- Mills equations are unsta- 
/\ • ble for arbitrary gauge groups. The proof involves the follow- 

ed ■ ing main steps. First, we show that the frequency spectrum of 

a class of radial perturbations is determined by a coupled sys- 
tem of radial "Schrodinger equations" . Eigenstates with negative 
eigenvalues correspond to exponentially growing modes. Using 
the variational principle for the ground state it is then proven 
that there always exist unstable modes (at least for "generic" 
solitons). This conclusion is reached without explicit knowledge 
of the possible equilibrium solutions. 
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1 Introduction 

In recent years the study of regular and black hole solutions of the Einstein- 
Yang-Mills (EYM) equations, without and with additional fields (Higgs fields, 
dilaton fields, etc.) has been actively pursued, and a number of interesting 
and surprising results have been discovered. Among the regular solutions 
the most interesting ones are those for which gravity is essential. The first 
example of this type was found by Bartnik and McKinnon [|l|] for the SU{2) - 
EYM system. For the same model several authors 0, ^ ^ discovered later 
the colored black hole solutions which showed that the classical uniqueness 
theorem for the Abelian case does not generalize. The existence of both types 
of solutions, which meanwhile has been established rigorously p, |], |^, led to 
a search for corresponding solutions of other related field theories. It turned 
out, for instance, that the Einstein- Skyrme system has black hole solutions 
with hair which are at least linearly stable P, ^, 0, |ll|. (For a numerical 



investigation of nonlinear stability, see ref. |TT1.) Several authors studied 
other models, notably the SU{2) -EYM system with a Higgs triplett |jl2|, |13, 



I^ , |T^, as well as the EYM-dilaton system [0, and found in some cases 



other linearly stable black hole solutions. Interesting black hole solutions 
have recently been found numerically for the EYM system with a Higgs 
doublet fl^, as in the standard electroweak model. These "sphaleron black 
holes" are expected to be unstable, but this question is not yet fully analyzed. 
We have shown recently, that the regular solutions are unstable |18|, but our 
method cannot be applied directly to the black hole case. 

The Bartnik-McKinnon solution and the related black hole solutions were 



shown to be unstable by some of us ||19|, ^, 0, ^. On physical grounds one 



would expect that this remains true for any gauge group, but a mathematical 
proof of this conjecture presents quite a challenge. We finally succeeded in 
constructing such a proof for the regular solutions, which will be sketched in 
the present paper. Similar arguments may perhaps also be used successfully 
to establish the instability of all black hole solutions, but problems related 
to the boundary conditions at the horizon have not yet been overcome. 

Our strategy is based on the study of the pulsation equations describing 
linear radial perturbations of the equilibrium solutions. Because the deriva- 



tion of these equations relies heavily on our previous work p3| , Q, we first 
recall some basic facts and equations in the next section. The pulsation 
equations are then presented in section 3 in a convenient partially decoupled 
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form. It turns out, that the frequency spectrum of a class of radial pertur- 
bations is determined by a coupled system of radial "Schrodinger equations" 
whose bound states correspond to exponentially growing modes. In section 
4 we prove with the variational principle for the ground state that there al- 
ways exist unstable modes, at least in the generic case (defined in sect. 2). 
The construction of the trial variations, which are used to establish this, is 
presumably the main point of the present letter; we succeeded only after a 
number of failures. 

2 Basic facts and equations 

Since we consider only spherically symmetric configurations, the metric g of 
the space-time manifold M can be parametrized as 

g = -NS^dt^ + N-^dr^ + r\d^^ + sin^^ dy^"^) , (1) 

where the metric functions A^ and S depend only on the Schwarzschild-like 
radius r and the time t. We use also the usual mass fraction m{r, t), defined 
by A^ =: l—2m/r. Gauge fields with spherical symmetry have been described 
in detail in ref. [^. Let us briefly recall the results, as far as they will be 
needed below. 

We fix a maximal torus T of G with the corresponding integral lattice 
I (= kernel of the exponential map restricted to the Lie algebra LT of the 
torus). In addition, we choose a basis S of the root system R of real roots, 
which defines the fundamental Weyl chamber 

K{S) = {H eLT\ a{H) > for all a G 5 } . (2) 

To a given spherically symmetric configuration, there belongs a canoni- 
cal element Hx G / n K{S) which characterizes the corresponding principal 
bundle P{M,G) admitting an SU{2) action (see sect. II of ref . P[). In ^ 



we showed that Hx is restricted to a small finite subset of / fl K{S), if the 
solution is regular at the origin. In much of our discussion we exclude the 
possibility that Hx lies on a boundary of the fundamental Weyl chamber and 
consider only what we call the "generic case", for which Hx is in the open 
Weyl chamber K{S). 



Instability of EYM Solitons 3 

A spherically symmetric gauge field can be described by a smooth family 
of linear maps A: LSU{2) — > LG, depending only on r and t and satisfying 
the following conditions: 

Ai = [A2, A3] , A2 = [A3, Ai] , A3 = -i/A/4vr , (3) 

where A^ := h{Tk) and 2iTk are the Pauli matrices. These equations imply 
that A+ := Ai + 1^2 hes in the following direct sum of root spaces La. 

aes (4) 

E:={aeR+\ a{Hy) = 2}. 

In the generic case E is a basis of a root system contained in R (see Appendix 
Aof ref. pi). 



For the example of the gauge group SU{2), H\ is an integer multiple of 
47r T^: H\ = Aiik t^ with k & Z^ and the only solutions of (^ are Ai = A2 = 0, 
A3 = fcrg, or 

Kl=WTi+WT2, A2 = TW Ti±WT2 , A3 = ±T3 . (5) 

In ref. [^ it is shown that there exists always a (local) gauge such that 
the gauge potential A takes the form 

A = A + A, (6) 

with 

i = A2rf^+ (A3Cos^- Aisin^)rfV5 (7) 

and 

A = NSAdt + Bdr, (8) 

where A and B commute with Hx (i.e. with A3). If Hx is generic one knows 
that its centralizer is the infinitesimal torus LT. Hence, A and B are LT- 
valued and A is thus Abelian. 

The coupled EYM equations, corresponding to the parametrizations (|l|), 
® ^ ® have been derived in pSf . Here, it suffices to write them in a slightly 
different form for the temporal gauge ^ = 0. 

The Einstein equations reduce to (k := SttG) 

m'='^{NG + pe}, m='^NH, (9) 
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where 



b r 



(10) 



G = -{{NS)-'\A+\'+ A'^ + [B,A+]\'\, 


(11) 


H = Re(A+, AV+[i3,A+]), 


(12) 


P^ = ^{-^'+-^11 


(13) 


^=^^, ^=^[A+,A_]-A3. 


(14) 



with 



Here we have used the following notation: ( ■ , ■ ) denotes a suitably nor- 
malized Ad(G)-invariant scalar product on LG, as well as its (hermetian) 
extension to LGc, and | ■ | is the corresponding norm. 
The YM equations decompose into 



^r^Hi + 



NS\S 



A+ , A'_ + [B, A_ 



A_,AV + [i3,A+] =0, (15) 



S\NS 



K) -^{NS{A'^ + [B,A^]}y 



- n{ [S, A'J + [B, [B,A^]]} + ^ [^, A+] = 



Ki ^) + 4f^' ^] + -ws^ ['^+'^-] + ['^-'^+] ) = ° • 



(16) 
(17) 



The last equation is the Gauss constraint. For the generic case the term 
proportional to [B, B ] in ([T7|) vanishes. 

For static solutions all time derivatives disappear and the basic equations 
simplify considerably. (For the Bartnik-McKinnon solution A is of the form 
(D with w = 0, A3 = rg and i = in (D.) 
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3 Pulsation equations 

We consider now a static, regular, asymptotically flat solution of the coupled 
EYM equations (|), (p!0D , (|15D - (|l^ with a generic Hx- Such a solution is of 
purely magnetic type with vanishing YM charge. This is not yet rigorously 
proven under satisfactory weak assumptions, but there is strong evidence 
for this (see |2^, ^ for partial results). From now on the symbols A±, N, 
S, etc. will be used to denote the equilibrium solution. Time-dependent 
perturbations are denoted by SA±, 6B, etc.. 

The details of the reduction of the perturbation equations to a well- 
adapted form will be described elsewhere. Here, we emphasize only some 
crucial points and present the final result. 

An observation which was already pointed out in [|l^ for G = SU{2) 
finds here a natural generalization which simplifies matters considerably: 
The first of the Einstein equations (P) leads to an expression for Sm' which 
can be readily integrated with respect to r. An arbitrary function of t which 
thereby appears for 6m is then fixed by the second Einstein equation in (|^). 
We find 

5m=-NRe{A'_^,SA+) . (18) 

It turns out that a significant decoupling of the perturbation equations 
is achieved by decomposing 6A^ as follows. Let us choose in the space (^ a 
basis {cq} of the root spaces La (a G S), with respect to which we expand 
the unperturbed A+ and its perturbation 5A+, 

6A+ = J2 Sw"ea. (19) 

Then we have ((5A_ := c((5A+), c = conjugation in LGc) 

6A± = 6X± ± i6Y± (20) 

with 

5X+ = ^ Re (5w") e„ , 5Y+ = Y. Im {Sw") e^ ■ (21) 

We shall call SX±, 6Y± the real and imaginary parts of the perturbations 
6A±. The unperturbed quantity A+ can be chosen to have only a real part, 
as was shown in [^ . 
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With the help of the equihbrium equations one can now bring the per- 
turbation equations after some work into the following standard form: 6X^ 
decouples, 

5X+ + Uxx SX+ = 0, (22) 

and 5$ := (^1+, SB) satisfies 

T5$ + U6^ = , (23) 

where the operators Uxx and U, 

t/=f^'"^ ^^"1 , (24) 

y UjSY Ubb j 

are given by the expressions 

TV 9^ 1 

Uxx = P*^ + -^^d{tJ^) - ^^ad(iVSA+)ad(iV^A_) 

+ (p*A+)^([-^,A+],-) + [^,A+]^(P*A+,-), (25) 

Uyy = P*^ + ^ad(i^) , (26) 

Ubb = - ad(iV5A+) ad(A^5A_) , (27) 

Uyb = +p,ad{NSA+) + iid{NSp,A+), (28) 

Uby = -ad(A^SA_)p, + ad(A^Sp,A_.) (29) 

with 

p, = -tNs4- ■ (30) 

ar 

The operator T in ( p3D acts as multiplication with the diagonal matrix 



T 



/ 



' "]. (31) 

Nr"^ J 
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It is easy to see that U and T are symmetric relative to the following 
scalar product for LGc-valued functions on 7?.+ : 

/■OO f]f 

(*i<^)=y^ (^,<^) ^. (32) 

(Questions of domains of unbounded operators will be discussed elsewhere.) 



For a harmonic time dependence, 5$ = ^e *'^*, we obtain from (p3D 



-^=^^- (33) 

{i\T\i} 

Different eigenmodes are orthogonal with respect to the following scalar prod- 
uct defined by T: (^i|^2)t := (^i| ^1^2)- For the frequency uq of the 
fundamental mode of eq. (|23|) we have the minimum principle 

^,2 ^ Min ill^jll , (34) 

which will be used in the next section to show that ujq^ < 0. 

The perturbation equations (^), (^) do not include the Gauss constraint 



(|T^, whose linearization reads 
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j^{p.[^SB^-[A^,6Y_]}=0. (35) 

It will turn out that this is automatically satisfied for physical pulsations 
(more precisely: in the space orthogonal to pure gauge modes). 

4 Instability of all generic EYM solitons 

For a given regular solution with A+ = J^aes W^^a we construct now a one- 
parameter family of field configurations A(j^)+, B(^y^) such that the variational 
expressions (p^ for 5K± = {dA(^^^^/dx)x=o^ ^^ = ('^'^(x)/'^x)x=o is negative 
. This family is chosen of the following form: 

A(^)+ = Ad(exp(-x^)){ A+ cos(x) + iA+(oo) sin(x) } , (36) 

%) = XZ', (37) 
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where Z is an LT-valued function of r with the properties 

hm [Z,A+]=iA+{oo) , suppZ' C [1 - e , 1 + e] (38) 

r— >0,oo 

for an e > 0. The existence of such a function can be seen as follows. Let 
{ha}a&T: be the dual basis of 27rS and put 

7_\-7"7, 7a _ / «^"(oo)/w"(0) forr<l-e, . . 

Z-l^ Z h^, ^ - 1 1 for r > 1 + e . ^"^^^ 

It is easy to verify that both conditions in (^) are satisfied. (It can be shown 
that m;"(0) ^ for all a G S; see Appendix A of ref. 0.) 

We note a few properties of the family (0), (^). The equilibrium solution 
is clearly obtained for x = 0. Applying a gauge transformation with g = 
exp{xZ) we obtain 

A(^)+ ^ A+ cos(x) + ^A+ (00) sin(x) , S(x) ^ . (40) 

The first variations of ( PB] ) and (^) are 

(5A+ = -[Z,A+]+zA+(oo), 6B = Z', (41) 

and these satisfy by construction the desired boundary conditions 

lim 5A+ = , lim 6B = . (42) 

r— >0,cx3 r^0,oo 

((5i3 vanishes even outside [ 1 — e , 1 + e ] .) Finally, (5A+ has only an imaginary 
component 

6Y+ = 5A+ = -[Z, A+] + iA+(oo) (43) 

and thus by (H) 

2 (6<^\U\6<^) 



with 5$ = (5y+,(5i3) given by (J3) and the second eq. in (41). 

It is easy to see that the denominator in this expression is finite if A_|_(r) 
A+(oo) vanishes sufficiently fast for r —>■ cxd. 



The calculation of the numerator in (|4^) for the operator U in ( P4|) is 
somewhat tedious. After several partial integrations, using ( ^2]) and the 
equilibrium equations, one finally arrives at 

{S^\U\S^) = -f {nS |A'+|' + 2^\:F\^ } dr , (45) 
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which is finite and strictly negative! 

We now return to the hnearized Gauss constraint (^) and note that a 
variation (5$ is orthogonal with respect to ( ■ | ■ )j, to all gauge variations 

'^^gauge = i^[x, A+] , x' ) , X-^^ ^ LT , (46) 

if and only if the curly bracket in (^) vanishes. This can be shown readily 
with a partial integration. 

This proves the instability of the regular equilibrium solution. It should 
be emphasized that we were able to draw this conclusion assuming only weak 
asymptotic conditions for the solitons. Especially, we have used repeatedly 
that the YM charge vanishes only if lim^^oo A(r) is a homomorphism from 
LSU{2) to LG (see eq. (|l|)). 
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